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Transient  surface  supersaturation  after  crystal  submersion  ( II ) 


Miroskwa  Rak* 

Institute  of  Physics,  Technical  University  of  Lodz,  ul.  Wolczanska  219,  93-005  Lodz,  Poland 


ABSTRACT 

In  our  earlier  papers  [ M.  Rak,  Proc.  SPIE  3 178  (1997)  108  and  M.  Rak  et  al.,  J.  Cryst.  Growth  197  (1999)  944  ],  we  solved 
the  time-dependent  equation  of  Burton,  Cabrera  and  Frank  (BCF)  for  transient  state  conditions  occurring  immediately  after 
submersion  of  a crystal  surface  in  supersaturated  solution.  As  a result,  the  expression  for  the  transient  surface 
supersaturation  was  found  and  discussed.  In  this  paper  we  present  an  approximation  of  that  expression.  The  approximation 
enables  a simple  estimation  of  the  time  required  to  attain  the  steady-state  value  of  the  surface  supersaturation. 

Keywords:  time-dependent  BCF  equation,  surface  supersaturation,  surface  diffusion  theory  of  crystal  growth,  solution 
growth 


1.  INTRODUCTION 

In  our  previous  papers  ’’ 2 we  studied,  using  the  surface  diffusion  theory  3’ 4 of  Burton,  Cabrera  and  Frank  (BCF),  the 
transient  surface  supersaturation  as  occurring  immediately  after  submersion  of  crystal  face  in  supersaturated  solution.  For 
this  purpose,  we  solved l>  2 the  BCF  time-dependent  equation: 
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where  n = n(x,t)  is  the  local  concentration  of  growth  units  at  the  surface;  n«  is  the  equilibrium  concentration  of  growth  units 
at  the  surface;  cr  is  the  relative  supersaturation  just  above  the  surface  and  very  far  from  a step,  so  it  is  the  same  as  the  bulk 
supersaturation;  X is  the  mean  diffusion  distance  of  the  growth  unit  adsorbed  on  the  surface  and  x denotes  the  relaxation 
time  for  leaving  the  surface  adsorption  layer.  The  relaxation  time  x is  related  to  the  activation  free  energy  AG  for  desorption 
by  Eyrings’  formalism3: 
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where  h is  Planck’s  constant,  kB  is  the  Boltzmann  constant  and  T denotes  temperature. 


To  formulate  the  initial  condition,  we  assumed  *’ 2 that,  as  soon  as  the  crystal  face  is  submerged  in  the  supersaturated 
solution,  the  surface  supersaturation  cts  is  equal  to  the  bulk  supersaturation  o: 


as(x,0)  = a , 


(3) 


where  the  surface  supersaturation  cts  ■ crs  (x,  t)  is  defined  as: 
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Eq.  (1)  was  solved  ’’ 2 for  a parallel  sequence  of  equidistant  steps  and  the  origin  of  the  coordinate  system  chosen  in  the 
centre  of  the  step  terrace.  The  positions  x = ± xo/2  (where  xo  denotes  the  interstep  distance)  corresponding  to  the  step  edges 
were  assumed  *■ 2 to  remain  unchanged  during  the  transient  state.  Furthermore,  the  boundary  conditions  assumed  T>  2 
symmetrical  diffusion  fields  around  steps  and  equilibrium  concentration  n of  growth  units  at  the  step  edges.  Taking  into 
account  Eq.  (4),  the  solution  for  Eq.  (1)  with  initial  condition  (3)  was  found 1,2  in  the  following  form: 
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The  analysis  of  the  transient  surface  supersaturation,  expressed  by  Eq.  (5),  revealed 1,2  that  the  time  ts  required  to  attain  the 
steady-state  value  of  cts  does  not  depends  on  the  bulk  supersaturation  a,  but  it  increases  with  increasing  relative  interstep 
distance  XcA.  Furthermore,  the  longest  period  of  time  tso  is  necessary  to  reach  the  steady-state  value  of  os  at  the  centre  of 
the  step  terrace  (x  = 0)  and  therefore  this  point  is  of  particular  interest 


The  aim  of  this  paper  is  to  present  an  approximation  of  Eq.  (5)  which  enables  a simple  estimation  of  the  time  ts  required  to 
attain  the  steady-state  value  of  surface  supersaturation  after  submersion  of  crystal  face  in  the  supersaturated  solution. 


2.  APPROXIMATION  OF  EQ.  ( 5 ) 

In  order  to  find  an  approximation  of  Eq.  (5)  we  follow  in  the  way  similar  to  that  presented  in  one  of  our  earlier  papers 5. 
Therefore,  we  use  the  effective  time  t*®  defined 5 as: 


cosh(x/X) 
cosh(x0  / 2X) 
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In  the  summation  term  of  Eq.  (5),  we  now  replace  the  power  f 7t2(X/x0)2(2k- 1 )2+ 1 jt/x  by  Wx,  which  leads  to  an 
approximation: 
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Taking  into  account  that  the  solution  given  by  Eq.  (5)  fulfils  the  initial  condition  [Eq.  (3)],  we  find: 
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This  result  is  also  obtainable  using  the  Fourier  expansion  corresponding  to  the  function  cosh(xA,)/cosh(xo/2A.)  in  the  interval 
(-x</2 , Xq/2).  Upon  substituting  Eq.  (8)  into  Eq.  (7),  we  have  an  approximation  of  Eq.  (5)  in  the  following  form: 
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where  W is  expressed  by  Eq.  (6).  The  approximation  given  by  Eq.  (9)  and  the  solution  found  previously  [ Eq.  (5)  ] take  the 
same  initial  and  steady-state  values.  In  addition  to  this,  Eq.  (9)  fulfils  the  boundary  conditions  assumed  for  solving  the  BCF 
time-dependent  equation. 


3.  DISCUSSION 

In  order  to  illustrate  the  solution  given  by  Eq.(5)  and  its  approximation  [Eq.  (9)],  we  perform  calculations  for  bulk 
supersaturations  sufficiently  low  to  predict  the  BCF  spiral  growth  mechanism  The  relative  interstep  distance  xoA.  present  in 
Eqs.  (5)  and  (6)  can  be  found  using  one  of  the  known  approximations  of  the  interstep  distance  (see  e.g.  Ref 6 ).  However, 
since  this  distance  depends  not  only  on  supersaturation  but  also  on  the  edge  free  energy  and  the  number  of  cooperating 
spirals,  we  assume  several  values  of  xq//-  (xo/k  = 2,  3,  4,  5)  to  show  the  influence  of  Xq/a  on  behaviour  of  the  surface 
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supersaturation.  In  Fig.  1 the  surface  supersaturation  cso  at  the  centre  of  step  terrace  (x  = 0)  is  shown  as  a function  of  the 
relative  time  t /x. 


t h 

Fig.  1.  Surface  supersaturation  aso  at  the  centre  of  the  step  terrace  (x  = 0)  versus 
relative  time  t/x  for  bulk  supersaturation  a = 0.01.  Curves  1,  2,  3 and  4 
correspond  to  the  relative  interstep  distance  x<A  = 2,  3,  4 and  5,  respectively. 

Full  and  dashed  lines  correspond  to  Eq.  (5)  and  (9),  respectively. 

The  full  and  dashed  lines  correspond  to  oso  calculated  from  Eq.  (5)  and  (recalculated  from  Eq.  (9),  respectively.  It  can  be 

seen  in  Fig.  1 that,  at  small  values  of  t/x,  the  function  aSo(t)  differs  from  the  exponential  function  c& , (t) . The  difference  is 

more  considerable  for  greater  values  of  the  relative  distance  xofk.  The  slope  (das/dt\ = 0 of  the  curve  as(t)  at  t = 0 is  equal  to 
zero,  which  means  that  at  t = 0,  the  surface  supersaturation  cs  remains  constant  (cf.  Ml  lines  in  Fig.  1).  This  result  is 
obvious  from  physical  point  of  view.  At  t = 0,  the  behaviour  of  the  surface  supersaturation  is  determined  by  the  exchange  of 
growth  units  between  the  crystal  surface  and  the  solution  bulk.  The  net  flux  jv  of  growth  units  exchanged  with  the  bulk  of 
the  solution  is  equal 3’ 4 to: 

jv  = ne(a-Qs)  (10) 

x 

However,  just  after  submersion  of  the  crystal  face  in  supersaturated  solution,  the  surface  supersaturation  crs  is  equal  to  bulk 
supersaturation  ct  [cf.  Eq.  (3)].  In  consequence,  at  t = 0,  the  net  flux  jv  is  equal  to  zero  and  crs  remains  constant  at  t = 0.  At 
t > 0,  the  surface  supersaturation  is  reduced  by  the  surface  flux  of  growth  units  going  to  the  step  edges. 

In  Fig.  2 the  surface  supersaturations  as  and  Oj  are  shown,  for  xq/a  = 3 and  various  values  of  t/x,  as  a function  of  the 
relative  distance  x/xo  from  the  step  centre.  It  can  be  noticed  that  with  increasing  distance  x/xq  from  the  step  centre,  the 
difference  between  cts  and  Cg  diminishes  and  for  x = (0.2  + 0.3)xo  reaches  zero.  Then  it  changes  its  sign  and  with 

increasing  x,  the  absolute  value  of  the  difference  increases.  Consequently,  for  x = (0.2  -s-  0.3)xo,  values  of  crs  and  c$  are 
very  close  to  each  other. 

Although  the  behaviour  of  the  functions  crs(t)  and  a*  (t)  is  different  at  the  beginning  of  the  transient  state,  these  both 
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Fig.  2.  Surface  supersaturation  cts  versus  relative  distance  x/xo  from  the  terrace  centre, 
plotted  for  XoA.  = 3,  <r  = 0.01  and  various  values  of  the  relative  time  t/x.  Curves 
1, 2,  3 and  4 correspond  to  t/x  = 0.1, 0.5,  1 and  3,  respectively.  Full  and  dashed 
lines  correspond  to  Eq.  (5)  and  (9),  respectively. 

functions  attain  their  steady-state  value  at  the  same  time  (cf.  Fig.  1).  Therefore,  the  time  ts  required  to  attain  the  steady-state 
value  of  the  surface  supersaturation  can  be  estimated  with  the  use  of  Eq.  (9).  In  transient  state  analysis,  it  is  usually  assumed 
that  the  steady-state  is  attained  when  the  transient  part  of  the  equation  is  not  greater  than  0.01  of  its  initial  value. 
Consequently,  the  exponential  term  in  Eq.  (9)  can  be  neglected  when  it  is  not  greater  than  0.01  of  its  value  at  t = 0.  This 
means  that  the  time  ts  necessary  to  attain  the  steady-state  value  of  the  surface  supersaturation  can  be  estimated  as  . 


ts«4.6t*ff  , (11) 

where  teg  is  given  by  Eq.  (6).  Values  of  the  characteristic  time  tso  (for  x = 0)  and  the  time  ts  for  x = 0.25xo,  obtained  using 
Eq.  (11),  are  presented  in  Table  1. 


Table  1.  Values  of  the  relative  time  ts/x  obtained  using  Eq.  (1 1)  and  those  following  from  Eq.  (5),  for  various 
values  of  xq/X.  Results  for  x/xo  = 0 and  x/xq  = 0.25. 


XoA. 

Wx 

Eq.(6) 

ts/x 

Eq.(ll) 

ts/x 

2 

0.35 

1.62 

1.43 

3 

0.57 

2.64 

2.41 

4 

0.73 

3.38 

3.22 

5 

0.84 

3.85 

3.87 

6 

0.90 

4.14 

4.41 

7 

0.94 

4.32 

4.89 

8 

0.96 

4.43 

5.32 

XoA. 

Wx 

Eq.(6) 

ts/x 

Eq.(ll) 

ts/x 

2 

0.27 

1.24 

1.29 

3 

0.45 

2.07 

2.12 

4 

0.59 

2.71 

2.74 

5 

0.69 

3.18 

3.17 

6 

0.77 

3.53 

3.47 

7 

0.82 

3.78 

3.70 

8 

0.86 

3.97 

3.88 
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In  this  table,  we  also  show  values  of  the  characteristic  time  tso  and  the  time  ts  calculated  from  Eq.  (5).  To  calculate  the  times 
tso  and  ts  from  Eq.  (5)  we  also  assumed  that  the  transient  part  (summation  term)  in  Eq.  (5)  can  be  neglected  when  it  is  not 
greater  than  0.01  of  its  initial  value.  It  follows  from  the  results  presented  in  Table  1 that,  in  the  centre  of  step  terrace  (x  = 0), 
the  simple  expression  (11)  may  be  used  for  estimation  of  the  characteristic  time  tso  only  at  not  too  large  values  of  the 
interstep  distance  (xoA.  < 7).  At  large  xq/X,  values  of  tso  estimated  using  Eq.  (11)  are  too  low  in  comparison  with  those 
following  from  Eq.  (5).  It  can  also  be  noted  that  at  the  point  x=0.25xo,  where  the  difference  between  cs  and  Cg  is  negligible 
(cf.  Fig.  2),  the  difference  between  values  of  ts  estimated  using  Eq.  (11)  and  those  following  from  Eq.(5)  is  also  negligible. 

4.  CONCLUSIONS 

The  present  study  shows  that  for  not  too  large  interstep  distances  XoA,,  the  time  ts  necessary  to  attain  the  steady-state  value 
of  the  surface  supersaturation  cts  is  approximately  equal  to  4.6  Ue,  where  W is  expressed  by  Eq.  (6).  This  simple  estimation 
demonstrates  that  the  time  ts  does  not  depends  on  value  of  bulk  supersaturation  but  it  increases  with  increasing  interstep 
distance  XoA.  and  with  increasing  distance  from  the  step  edge.  It  should  also  be  pointed  out  that  the  time  ts  is  directly 
proportional  to  the  relaxation  time  x for  leaving  the  surface  adsorption  layer.  This  means  that  the  time  ts>  required  to 
attain  the  steady-state  value  of  os,  is  related  to  the  activation  free  energy  AG  for  desorption,  and  temperature  of  the 
solution  [cf.  Eq.  (2)]. 


5.  LIST  OF  SYMBOLS 

AG  Activation  free  energy  for  desorption, 

jv  Net  flux  of  growth  units  exchanged  with  the  solution  bulk, 
h Planck’s  constant. 

kB  Boltzmann  constant 

n [=  n(x,t)]  Local  concentration  of  growth  units  at  the  surface, 

n,  Equilibrium  concentration  of  growth  units  at  the  surface. 

T Temperature, 

t Time. 

Lff  Effective  time  defined  by  Eq.  (6). 

ts  Time  required  to  attain  the  steady-state  value  of  os. 

tso  Characteristic  time,  required  to  attain  the  steady-state  value  of  the  surface  supersaturation  Oso- 

x Distance  from  the  step  centre. 

Xo  Interstep  distance. 

X Mean  diffusion  distance  of  the  growth  unit  adsorbed  on  the  surface. 

a Relative  supersaturation  just  above  the  surface  and  very  far  from  a step,  so  it  is  the  same  as  the  bulk 

supersaturation 

cts  [-  as(x,t)]  Surface  supersaturation  expressed  by  Eq.  (5) 

Oj  [s  a*  (x,  t)  ] Surface  supersaturation  cts  approximated  by  Eq.  (9). 

<Tso  Surface  supersaturation  cts  at  the  centre  of  the  step  terrace  (x=0). 
a g0  Surface  supersaturation  ctso  approximated  by  Eq.  (9). 
x Relaxation  time  for  leaving  the  surface  adsorption  layer. 
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